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We numerically calculate the nuclear spin-lattice relaxation rate (Rs) in the superconducting 
state in terms of anisotropic s-wave gaps. By taking into account electron-phonon coupling, our 
calculated R s for a conventional s-wave superconductor, indium, is in quantitative agreement with 
the experimental data with a clear Hebel-Slichter peak. In contrast, by using the highly anisotropic 
s-wave gaps inferred from the magnetic penetration depth and scanning tunneling microscopy, our 
calculated R s curves for electron-doped Pro.g1LaCeo.09CuO.4_,, show no Hebel-Slichter peak, in 
agreement with the experimental data. Finally, the observed weak nonmagnetic pair-breaking effect 
provides unambiguous evidence for a highly anisotropic s-wave gap in this underdoped cuprate. 



The identification of the intrinsic gap symmetry in 
cuprates is crucial to the understanding of the micro- 
scopic pairing mechanism of high-temperature supercon- 
ductivity, which remains elusive for over twenty years. 
The superconducting transition temperatures T c 's of 
hole-doped cuprates appear to be too high to be ex- 
plained by the conventional phonon-mediated pairing 
mechanism. In contrast, the highest T c in electron-doped 
(n-type) cuprates is about 40 K, which is within the 
T c limit of the conventional phonon-mediated mecha- 
nism. Indeed, earlier [lj and recent Q tunneling spectra 
in electron-doped cuprates show strong electron-phonon 
coupling features, similar to the conventional supercon- 
ductors. The predominantly phonon-mediated pairing 
should be compatible with an s-wave gap. Many inde- 
pendent experiments designed to test the gap symme- 
try in the electron-doped system have led to controver- 
sial conclusions. Surface-sensitive angle-resolved photoe- 
mission spectroscopy (ARPES) [!, |j| implies a d-wave 
gap with a maximum gap size of about 2.5 meV. This 
gap size would imply a T c of about 14 K at the top 
surface, which is a factor of 1.9 lower than the bulk 
T c of 26 K (Ref. [4J). Surface and phase-sensitive ex- 
periments Q provide evidence for pure d-wave order- 
parameter (OP) symmetry in optimally doped and over- 
doped n-type cuprates. In contrast, nearly bulk-sensitive 
point-contact tunneling spectra along the Cu02 planes 
Q show no zero-bias conductance peak (ZBCP) in op- 
timally doped and overdoped samples [7Hllj. which ar- 
gues against d-wave gap symmetry. The bulk-sensitive 
Raman scattering data of Ndi.ssCeo.isCuO^y imply an 
anisotropic s-wave gap with a minimum gap A m i n of 3.2 
meV (Ref. [13): which is very close to A m in — 3 meV 
inferred from the magnetic penetration depth data [Hj]. 
Bulk-sensitive thermal conductivity [l4| and specific heat 
p~5j | data seem to support d-wave gap symmetry [fl HH 
while the same data can be quantitatively explained by 
nodeless s-wave gap symmetry [Hi ]. The absence of the 
"Hebel-Slichter" or "coherence" peak below T c in the nu- 
clear spin-lattice relaxation rate (R s ) of a slightly under- 
doped n-type Pr .giLaCeo.o9Cu0 4 - a (T c = 24 K) (l7j 



appears to argue against s-wave gap symmetry. How- 
ever, the absence of the coherence peak does not nec- 
essarily rules out s-wave gap symmetry because there 
are several mechanisms that can suppress the coherence 
peak. One of the mechanisms is quasi-particle damp- 
ing due to strong electron-phonon coupling [3]. Strong 
electron-electron correlation also leads to a strong sup- 
pression of the coherence peak [19]]. Furthermore, the 
coherence peak can be also reduced by gap anisotropy 
20l 21 1 . Therefore, the absence of the coherence peak 
in the R s data of Pro.giLaCeo.osjCuO^j, may arise from 
the combination of the intermediate electron-phonon cou- 
pling constant (A ~ 1) [lj], strong electron-electron cor- 
relation, and a highly anisotropic s-wave gap. 

Here we present numerical calculations of the nu- 
clear spin-lattice relaxation rate in the superconduct- 
ing state in terms of anisotropic s-wave gaps. By tak- 
ing into account electron-phonon coupling, our calcu- 
lated R s for a conventional s-wave superconductor, in- 
dium, is in quantitative agreement with the experi- 
mental data with a clear Hebel-Slichter peak. In con- 
trast, by using the highly anisotropic s-wave gaps in- 
ferred from the magnetic penetration depth and scan- 
ning tunneling microscopy, our calculated R s curves for 
Pro.9iLaCeo.ogCu04_ a show no Hebel-Slichter peak, in 
agreement with the experimental data. Finally, the ob- 
served weak nonmagnetic pair-breaking effect provides 
unambiguous evidence for a highly anisotropic s-wave 
gap in this underdoped n-type cuprate. 

The expression for the ratio of R s /R n of an anisotropic 
superconductor with a complex gap function is given by 
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where kg is the Boltzmann constant, f{uj) is the 
Fermi-Dirac distribution function, < N(T, uS) > and 
< M(T,lu) > are the respective Fermi-surface averages 
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of N(T,u>) and M(T, w), which are given by 



N(T,u) = Re{ 
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and 



M(T,w) = i?e[ 



A(r,n) j 

cj 2 - A 2 (T, fi) 
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where fi is the direction vector on the Fermi-surface. 
Taking into account electron-phonon coupling, N(T,uj) 
and M(T,v) can be approximated by [21 1 
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M(T,uj) = Re[ 
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where Ai(T,fi) = A (T)A(fi) is the real part of the 
gap function, A(f2) determines the gap anisotropy, and 
S is given by [2l| 



mrX . T 



l [r(n + -)C(n- 



3., fc B T 
2^2A (T) 



3/2 



^,2A (T) )3/2cxp( _ Ao(T)/fcsT)L (6) 



Here T and C are the T and Riemann £ functions, 
respectively. At low temperatures, the electron-phonon 
spectral function a 2 (ui)F(ui) varies as u) n . In the follow- 
ing numerical calculations for i? s , we will adopt "jellium" 
model where n = 2 (Ref. [Hj|). 

We first present the results of numerical calculations of 
R s for a conventional superconductor, indium. In the cal- 
culations, we take the realistic parameters for indium: A 
= 0.8URef. 0), A (T) = 1.9T c tanh[1.81(T c /T- 1) 1 / 2 ] 
(Ref. [13]), and 6 D = 111 K (Ref. 0). We use three sim- 
ple gap functions to characterize gap anisotropy: A(U) 
= 1 for isotropic s-wave gap (s-wave-1); A(fi) = (1 — 
0.02 cos 40) for a slightly anisotropic s-wave gap (s-wave- 
2); A(f2) = (1 — 0.9 cos 4</>) for a highly anisotropic s-wave 
gap (s-wave-3). Fig. 1 shows numerically calculated re- 
sults of R s for the three anisotropic s-wave gaps together 
with the measured R s for pure indium. It is remarkable 
that the calculated curve for the isotopic s-wave gap (s- 
wave-1) is slightly off from the experimental data while 
the curve for the slightly anisotropic s-wave gap (s-wave- 
2) almost coincides with the data. In contrast, when 
the gap becomes highly anisotropic (e.g., s-wave-3), the 
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FIG. 1: Nuclear magnetic spin-lattice relaxation rate {Rs) 
of pure indium (solid circls). The data are taken from 
Ref. pH| . The numerically calculated R B for three gap func- 
tions: A(fi) = 1 for isotropic s-wave gap (s-wave-1); A(fi) = 
(1 — 0.02 cos 4(j>) for a slightly anisotropic s-wave gap (s-wave- 
2); A(n) = (1 — 0.9 cos 40) for a highly anisotropic s-wave 
gap (s-wave-3). 



Hebel-Slichter peak is completely removed. This suggests 
that a highly anisotropic s-wave gap can completely sup- 
press the coherenec peak. The numerical results for in- 
dium thus provide important insight into how sensitively 
the coherence peak changes with the gap anisotropy. 

In order to address whether the absence of the co- 
herence peak in the R s data of Pro.giLaCeo.o9Cu04_ y 
can be also explained in terms of a highly anisotropic 
s-wave gap, we need to extract the gap functions from 
independent experimental results such as tunneling spec- 
tra and the in-plane magnetic penetration depth A a {,(T). 
Fig. 2a shows normalized tunneling conductance of a 
electron-doped Pr .88LaCeo.i2Cu04_ y crystal (T c = 24 
K). The normalized tunneling spectrum is reproduced 
from Ref. Q • The spectrum was taken on the top CuC>2 
plane using scanning tunneling microscopy (STM) (24j. 
We can numerically calculate the tunneling conductance 
using the following equation [25j |: 
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where 9 is the angle measured from the Cu-0 bonding 
direction, T is the life-time broadening parameter of an 
electron, and p(9—9o) is the angle dependence of the tun- 
neling probability and equal to cxp[— /3sin 2 (# — O )]. The 
solid line is the numerically calculated curve using T — 
0.70 meV, j3 — 5.3, 9q = tt/4, and an anisotropic s-wave 
gap function: Ai(0) = 4.9(|1.43cos20-O.43cos60|+O.2) 
meV. The finite life-time broadening parameter T of 
an electron may be caused by disorder and inhomo- 
geneities. This life-time broadening effect at zero tem- 
perature should be also incorporated into the above ex- 
pressions for R s by replacing ui with to — iT. 
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FIG. 2: a) Normalized tunneling conductance of a 
electron-doped Pro.8sLaCeo.i2Cu04-j/ crystal (T c — 24 K). 
The solid line is numerically calculated curve using an 
anisotropic s-wave gap function: Ai(0) = 4.9(| 1.43 cos 26 — 
0.43 cos 69\ +0.2) meV. b) Temperature dependence of l/X^b 
for Lai.9i3Ceo.o87Cu04_ !/ (T c = 28.5 K). The solid line is the 
numerically calculated curve using an anisotropic s-wave gap 
function: Ai(0) = 3.5(| 1.43 cos 26 - 0.43 cos 60| + 0.35) meV. 

Figure 2b shows the temperature dependence of 1/ A^ b 
for Lai.9i 3 Ceo.o87Cu04_ 2/ (T c — 28.5 K), which has a 
similar doping level as PrngiLaCeo.o9Cu04_ i ,. The data 
are digitized from Ref. 26]. We can numerically cal- 
culate the temperature dependence of A 2 b (0)/A 2 & (T) for 
an anisotropic gap function using the following equation: 
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Here E = y/e 2 + A? (9, T) . The solid line is the 
numerically calculated curve using an anisotropic s- 
wave gap function: Ai(T,0) = 3.5tanh[1.81(T c /T - 
l) 1 / 2 ](|1.43cos20 - 0.43cos66»| + 0.35) meV. If we 
assume that the superconducting gap is propor- 
tional to T c , then the bulk superconducting gap of 
Pro.91LaCeo.09CuO4._j, inferred from the penertation 
depth data is 2.95(|1.43 cos 29 -0.43 cos 60| +0.35) meV. 



Now we can use the inferred gap functions from both 
STM and magnetic penetration depth to calculate R s 
for Pro.9iLaCeo.o9Cu04_ a . In the calculations, we take 
the realistic parameters for the electron-doped cuprate: 
A = 1.0 (Ref. and 9 D = 384 K (Ref. 0). Since the 
R s data were taken in a magnetic field of 6.2 T which 
suppresses T c to about 20 K, we scale down the gap size 
proportional to T c . For the gap function extracted from 
the tunneling spectrum, the same life-time broadening 
parameter (r = 0.7 meV) in the tunneling spectrum is 
used to calculate R s . This parameter only influences R s 
at low temperatures and has little effect on R s close to 
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FIG. 3: Nuclear magnetic spin-lattice relaxation rate (Rs) in 
Pro.9iLaCeo.o9Cu04_ !/ (solid circles). The solid and dashed 
lines are the numerically calculated R s curves for the s-wave 
gap functions inferred from STM (denoted as STM) and from 
the magnetic penetration depth (denoted as bulk). 

The calculated results are shown in Fig. 3. It is appar- 
ent that for the gap function extracted from STM, the 
calculated R s curve (solid line) agrees excellently with 
the experimental data. On the other hand, for the bulk 
gap function, the calculated curve (dashed line) deviates 
significantly from the data. In all the cases, the coher- 
ence peaks are absent. Therefore, the absence of the 
coherence peak in the R s data could be consistent with 
highly anisotropic s-wave gaps. 

Although the R s data can be excellently explained 
in terms of the gap function inferred from STM, it 
does not necessarily imply that the gap function de- 
duced from STM is more realistic than that inferred 
from the magnetic penetration depth. As a matter of 
fact, the above calculations do not take into account 
strong electron-electron correlation. It was shown that 
strong electron-electron correlation can reduce the co- 
herence peak of R s for an s-wave gap and cause R s to 
drop more rapidly just below T c j 1 01 ] - For the slightly 
underdoped Pro.9iLaCeo.ogCu04_ i ,, the measured value 
of T\TK 2 (where T\ is the spin-lattice relaxation time 
and K s is the spin part of Knight shift) is a factor of 50 
smaller than the expected value of noninteracting elec- 
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trons [Tjj. This implies strong electron-electron correla- 
tion. Therefore, it is very likely that the R s data can be 
quantitatively explained in terms of the bulk anisotropic 
s-wave gap if the strong electron-electron correlation is 
taken into account. 
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FIG. 4: The calculated curve of T c versus residual resis- 
tivity in terms of any d-wave gap (dotted line) and the 
bulk anisotropic gap (solid line): A = 2.95( 1 1.43 cos 26 — 
0.43 cos 6t9 + 0.35) meV. In the calculations, we use 
hflp = 0.62 eV. The solid circle is a data point for 
Pr .9iL a Ceo.o9Cu04- s/ (PLCCO) [13]. 



In order to unambiguously distinguish between any d- 
wave and anisotropic s-wave gap symmetries, we study 
the response of a superconductor to nonmagnetic impuri- 
ties or disorder. The nonmagnetic impurity pair-breaking 
effect is both bulk- and phase-sensitive. This is because 
the rate of T c suppression by nonmagnetic impurities [28j 
is determined by the value of the Fermi surface (FS) aver- 
age < A(fc) >fs, which depends sensitively on the phase 
of the gap function. More specifically, the rate is pro- 
portional to a parameter % = 1 — (< A(fc) >fs) 2 / < 
A 2 (fc) >fs- It is easy to show that \ — 1 for an Y <i-wave 
gap and \ = 0.058 for the highly anisotropic s-wave gap: 
A = 2.95(|1.43 cos 2(9-0.43 cos 66»| +0.35) meV. An equa- 
tion to describe the pair-breaking effect by nonmagnetic 
impurities (or defects) is given by [28| 
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where hfl* is the renormalized plasma energy [2 a. |29| in 
units of eV, p r is the residual resistivity in units of pQcm, 
and ^ is the digamma function. The lower limit of hfl* 
is equal to h£l* which is related to the zero-temperature 
penetration depth A a {,(0). For Lai.9i3Ceo.o87Cu04_j,, 
A q6 (0) = 320 nm (Ref. leading to hn* s = 0.62 eV. 

Figure 4 shows the calculated curves of T c versus resid- 
ual resistivity in terms of any d-wave gap (dotted line) 
and the bulk anisotropic s-wave gap (solid line): A = 
2.95(| 1.43 cos 2(9 - 0.43 cos 6(9 1 + 0.35) meV. For the d- 
wave gap, T c is suppressed to at p r — 74.6 ptlcm 



Pri.85Ce .i 5 CuO 4 - y (T c = 20 K), hti* is found to be 
a factor of 1.65 larger than HQ* due to a finite mean- free 
path [30l |. If we increase hfl* of Pro.9iLaCeo.o9Cu04_ y 
by the same factor (1.65), T c will be suppressed to 
at p r = 27 'A pVLcm for the d-wave gap. The measured 
large p r of 92 /zficm and nearly optimal T c of 24 K in 
Pro.9iLaCeo.o9Cu04_ i , rules out any <i-wave gap symme- 
try. 

In summary, the absence of the Hebel-Slichter 
peak in the nuclear spin-lattice relaxation rate of 
Pro.giLaCeo.o9Cu04_ a can be well explained by a highly 
anisotropic s-wave gaps inferred respectively from the 
magnetic penetration depth and scanning tunneling 
microscopy. The observed weak nonmagnetic pair- 
breaking effect provides unambiguous evidence for a 
highly anisotropic s-wave gap in this slightly underdoped 
n-type cuprate. 
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